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Abstract 

Supersymmetric quantum mechanics is a powerful tool for generating exactly solv¬ 
able potentials departing from a given initial one. In this article the first- and second- 
order supersymmetric transformations will be used to obtain new exactly solvable po¬ 
tentials departing from the complex oscillator. The corresponding Hamiltonians turn 
out to be ruled by polynomial Heisenberg algebras. By applying a mechanism to reduce 
to second the order of these algebras, the connection with the Painleve IV equation is 
achieved, thus giving place to new solutions for the Painleve IV equation. 
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1 Introduction 

As is well known, supersymmetric quantum mechanics (SUSY QM) is one of simplest tech¬ 
niques for generating new exactly solvable potentials departing from a given initial one [l]-|l0|. 
In general, the spectrum of the built Hamiltonian differs from the initial one in a finite num¬ 
ber of levels. In the past this procedure has been applied mainly to hermitian Hamiltonians, 


in order to generate either real or complex new potentials 111T , with an associated spec 


trum which in the last case could contain a finite number of complex energy levels additional 
to the initial eigenvalues. The method has been used also for non-hermitian Hamiltonians 
with the so-called PT-symmetry 


18 


Indeed, in the context of optics it is possible to 
find families of complex potentials with entirely real spectra with a missing single level, the 
one associated to an arbitrary initial mode 21 . However, up to our knowledge there are not 


many works employing Hamiltonians with purely complex eigenvalues. 


Our aim in this work is to address this subject. In order to do that, we will take as initial 

which has the form of the harmonic oscillator 


system the so-called complex oscillator 22-25 


potential but with a frequency u> which is allowed now to be complex. We are going to show 
that some particularly interesting first-order SUSY transformations, which in the real case 
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would introduce singularities in the generated potential, here will be non-singular and they 
will generate gaps in the equidistant spectrum of the complex oscillator. 


Let us recall that the SUSY partners of the harmonic oscillator potential supply explicit 
realizations of the even-order polynomial Heisenberg algebras (PHA) [26-33 


Moreover, 


there are subsets of this general family of SUSY partner potentials which are ruled simul¬ 
taneously by second-order PHA. As the last are directly linked with the Painleve IV (PIV) 
equation, then a simple technique is available for generating solutions to this equation de¬ 
parting from the extremal states of the Hamiltonian 31-40 . In fact, plenty of real as well 


32- 

38 

up 

39[ 


In 

all 


as complex solutions associated with real PIV parameters have been generated 
addition, several complex solutions for complex PIV parameters have been built up 
of them departing from the harmonic oscillator potential, which is real (for a review see, 
e.g. |401). However, since the complex oscillator can be mapped into the harmonic oscillator, 
thus we can generate complex solutions to the PIV equation directly by applying first-order 
supersymmetric transformations to the complex oscillator, and some special higher-order 
ones. 


This paper is organized as follows. In section [2] we will sketch the SUSY QM. Then, 
in section [3] we will discuss the PHA, with special emphasis placed in those of second- 
order, which have a direct link with the PIV equation. In section [4] we will generate the 
SUSY partners of the complex oscillator, while in section [5] we will study the corresponding 
solutions to the PIV equation. The final section contains our conclusions and some discussion 
concerning the perspectives of our work. 


2 Supersymmetric quantum mechanics 


Let us start with two Schrodinger type Hamiltonians 

1 d 2 


Hj 2 da; 2 


+ Vj (x), j = 0,1, 


and suppose the existence of two first-order differential operators A\ 




which intertwine Ho with the new Hamiltonian H\ in the way [6,7,41,42 

HiAt = A+H 0 , H 0 A^ = A^ Hi, 


( 1 ) 


( 2 ) 


(3) 


where f3\ (x, e) is an unknown complex function, so that Af and A 1 are not adjoint to each 
other. 

From Eqs. 0-6. the following relations have to be satisfied 

P'i ( x , e ) + Pi Ou e) = 2 [V 0 (x) - e ], 

V\ (x) = V 0 (a;) - P'i [x, e ), 
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(4) 

(5) 























with the same terminology as in the real case being used 126,43,44], i.e., e G C is the so-called 
factorization energy. Through the substitution 

A ( x , e) = [In u (x, e)]' = - f X ’ £ I , (6) 

U[X,£) 

Eq. Q is transformed into 

H 0 u(x,e) = eu (x, e) , (7) 

namely, u (x, e) is a solution of the stationary Schrodinger equation for H 0 associated with 
£. 


It is interesting to observe that the involved Hamiltonians become factorized as 

Hq = A.-^ Af + £, Hi = A^ Ai + £. ( 8 ) 

From the intertwining relationships (3), we can obtain the eigenfunctions 'ipn'* of Hi and 
their corresponding eigenvalues, departing from the normalized eigenfunctions '0n°' ) of Hq 
with eigenvalues E n : 


{x) = CnA+'l/jW (x) , 

Hi^ (x) = E n ipH (x), n = 0,1, 2, • ■ ■ , 


(9a) 

(9b) 


which means that 'ijjn' 1 and share the eigenvalues E n . The normalization factors C n in 
Eq. (9a) cannot be straightforwardly obtained, as for real u, since now A+ ^ (Aj")f. 


Note that there exists a wavefunction annihilated by A 1 given by 

1 


(®) 


oc 


u (x, £) ’ 


( 10 ) 


which satisfies Hi = e ■ If can b e normalized, it becomes also an eigenfunction 
of H i with eigenvalue e; in such a case, the spectrum of Hi is the same as for Hq plus an 
extra energy level at e: 

S P {Hi) = Sp (Hq) U {e}. 

In order to avoid singularities in the new potential Vi, the transformation function u ( x , e) 
(also called seed solution), must not have zeros. This is simple to achieve due to the com¬ 
plex nature of the function itself: hardly the real and imaginary parts of u (x, e) will vanish 
simultaneously at the same point. In addition, in this case the oscillation theorem does not 
hold anymore 45 . Due to this, from now on we are going to suppose that u (x, e) has not 
zeros on the real axis. 


By iterating now k times the first-order transformation just discussed, a k-th order com¬ 
plex SUSY transformation is generated, where k new energy levels can be created so that 
the spectrum of the corresponding Hamiltonian becomes [7,[43| 44,46,47 

Sp (H k ) = {ei, • • • ,£ k , E n ; n = 0,1, • • • }. (11) 


3 








The standard SUSY algebra with two generators Q i, Q 2 introduced by Witten [48| (see 
also [49]) 


( 12 ) 


\QjiH S s\ QjH ss H ss Qj 0 , 

{Qj 1 Qk} QjQk QkQj ^jkH ss , k 1? 2, 

can now be realized in the way 

Qi = 4 ( 0 + + Q -), Q2 = Y (Q+ _ Q -), 


where 


c + = 

^ = {Q",Q + } = 


/0 B+' 


1 ^ 

Vo o y 

), Q~ 

'B+B k 

0 

v 0 

B k Bt 


(H k 

H d = 



0 0 
0 


+ — (H d — £1) ■ ■ ■ ( H d — e k ) , 


(13) 


The so-called supersymmetric Hamiltonian H ss is a polynomial of degree k-th in H dl since 
H 0 and H k are intertwined as 


with 


being k-th order differential intertwining operators. The initial and final potentials, Vo and 
14 respectively, are related by 


H k B + = 

B+H„, H 0 B k = B k H t , 

(14) 

Bt = K 

■ ■ ■ Af , B k — A 1 ■ ■ ■ A k , 

(15) 


14 (z) = Vo (x) - [In W (ui, • • • , u k )]" , 

Uj = uj ( x) =: u (x, £j ), j = !,-■■ ,k, 


(16) 

(17) 


where VU (iti, • • • , u k ) is the Wronskian of the k Schrodinger seed solutions {wi, • • • , u k }. The 
eigenfunctions of H k associated to the eigenvalues E n are given by: 


(®) = C'n5+V’i 0) (z) , n = 0, 1,2, ■ ■ ■ 

On the other hand, those associated to the new eigenvalues e 3 become: 

W k ( Uj ) 


^ (x) OC , 

3 W (Ml, • • • ,u k ) 


3 = 1,’’’ , k , 


(18) 


(19) 


where W k {u 3 ) denotes the Wronskian of the seed solutions {ui, • • • , u k } excluding the jth one 
(uj). In particular, for k = 1 the first order case discussed previously is recovered if we assume 
that VU («i) = and W\ (iq) = 1 (compare Eqs. (5),(9a),(10) with Eqs. 

Let us note that this technique has been employed successfully to generate new solvable 
potentials V k departing from a given initial one V d for several interesting physical Hamilto¬ 
nians (mainly Hermitian 


26,43,44 
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3 Polynomial Heisenberg algebras 

A (m — l)-th order PHA is a deformation of the Heisenberg-Weyl algebra with three gener¬ 
ators {H, £+, } which obey two standard commutation relations 

[//,£*]= ±£±, (20a) 

and an atypical commutator characterizing the deformation: 

[£-, £+] =N m (H+ 1) - N m (H) = P m - 1 (if), (20b) 

where are m-th order differential ladder operators, N m (. H ) = £+£“ is a polynomial in 
H of degree m representing a generalization of the number operator, so that P m - 1 (if) is a 
polynomial in PI of degree m — 1. The generalized number operator can be factorized as 

m 

(H) = n ( H - A). (21) 

J=1 

Sj being the zeros of the polynomial N m (if), which correspond to the associated extremal 
state energies. 


Let us consider the kernel K of C m . i.e., the m-dimensional space of solutions of the 
m-th order differential equation 

C * = 0, (22) 


which is invariant under H. Then, from Eqs. (21) and (22) we obtain 


=n (h -=°- 

3 = 1 


(23) 


We choose as a natural basis of K the solutions of the stationary Schrodinger equation for 
H associated to £j G C: 

m e .=£ji) e .. (24) 

Departing from these so-called extremal states, by acting iteratively with £+ onto them, 
we can build m mathematical energy ladders with the same spacing A E = 1 starting from 
each £j. If s of these states are normalizable, {l/A, ; j — 1, • • • , s}, then Sp ( H ) will have s 
independent endless physical ladders (see Fig. [jji). It could happen, however, that for some 
j G {1, • • • , s} 

(At)”" 1 fe, + 0, (£+)”*, =0, n e N. (25) 

Then, by analyzing the expression £“(£,+ ) n ^ i = 0, we observe that one of the remaining 
non-physical roots must satisfy £ ( = £j + n, l E {s + 1, • • • , m}. Hence, Sp (H) will have 
s — 1 endless ladders and a finite one of length n, which starts from £j and ends at £j + n — 1 
(see Fig. 0>) • 

It is important to identify now the most general one-dimensional Schrodinger Hamilto¬ 
nians which are ruled by the polynomial Heisenberg algebras of Eqs. (18) for growing values 


5 




• £i + l 

• fi 


£ £2 + 1 

• £2 

(a) 


• £, + 1 

# £ s 


• fi + i 

• £1 


• Sj + n-l 

% £’2 + 1 

• £j+ 1 

• £2 



• £ s + l 


• £. 


Figure 1: Spectrum of the Hamiltonian H having either: (a) s physical extremal states from 
which s independent endless ladders are generated, or (b) s physical extremal states, s — 1 
of which are associated to endless ladders and a finite one with n steps. 


in the order m of the ladder operators Cf n 33 . The answer runs as follows: for m — 1 and 


m — 2 (ladder operators of first- and second- orders respectively) the most general potentials 
have a closed form, becoming the harmonic and effective radial oscillators respectively. On 
the other hand, for m = 3 (third-order ladder operators) there is not a closed expression for 
the potential, but it depends on a function g which must satisfy the PIV equation. Because 
its importance for this work, let us review quickly this last case. 

3.1 Second order PHA. 


The connection between PIV equation and PHA appears for rn — 3 33,49 . Thus, let 


£3 be a third-order differential ladder operator factorized as the following product of 
a first- and a second-order differential one: 


r+ — /’+/’+ 

- t -3 — -‘-'l z -2 5 


£0 = - 


d 2 d 

ds- +gix) ± c + h[x) 


(26a) 

(26b) 

(26c) 


/, g and h being complex functions to be determined. We assume the existence of an 
auxiliary Hamiltonian H a which is intertwined with H in the way 


HC+ = £+ (H. d + 1), H a C+ = C+H. 

After some work, the resulting uncoupled system of equations becomes 


(27) 


/ 

h 

V 


x + g, 



— 2V + £\ + £2 — 2, 


x 2 

T 



+ xg + £3 — 


1 

2 ’ 


(28a) 

(28b) 

(28c) 
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where g satisfies 


0" = |~ + ^9 3 + 4 xg 2 + 2 (x 2 - a) g + 
a := £i + So — 2£ 3 — 1 , b := —2 (£d — £ 2 ) , 


(29a) 

(29b) 


which is the PIV equation with parameters a,b G C [50] . Note that it is essential 
to know the function g associated to a pair of parameters a, b ; then, the extremal 


energies Sj, which now are complex, as well as the potential V in Eq. (28c) can be 


found. Moreover, the three extremal states become expressed in terms of g as: 


2 



flh+t-w- 


dx 


J ( 2 g + 2 + 9 \J ~2 ) d:r 


(30a) 

(30b) 

(30c) 


In addition, the generalized number operator is cubic in the Hamiltonian: 

N 3 (H) = (H- Si) (H - S 2 ) (H - S 3 ). (31) 

Conversely, if we are able to identify Hamiltonians having third-order differential ladder 
operators, it is possible to design a simple mechanism to obtain solutions to the PIV 
equation (see [33] and section [5] in this paper). 


4 SUSY partners of the complex oscillator 

In order to implement the SUSY technique, we need to solve first the stationary Schrodinger 
equation for the complex oscillator potential: 


V„(x)=iu 2 x 2 , ui = e' e , 


(32) 


with an arbitrary complex factorization energy e, u being a dimensionless complex frequency 


and 9 its phase. A direct calculation leads to 51 
u (x, e) = e _ 2 w 


1 9 

— e 2 UX 


4 2w'T 
1 £ 1 


1 A) ( 7 - ) + A Xi-Fi ( ^ ux 1 


1-^1 ( 7 + —, o! — UJX ) + A x iFi ( - + —, —< ujx‘ 


where we will take 


4 2u 2 

r 

1 V 4 2oj / 


3 _ £ 3 

4 ~ 2c? 2’ 
3 e3 


(33) 


X = 2u 


1 V 4 2oj / 


ll/l < 1, 
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and i F\ is the confluent hypergeometric function (Kummer). In general, u (x, e) is a divergent 
solution of the Schrodinger equation since it does not satisfy limi^i^oo | u(x,e) | = 0. To 
ensure its square-integrability, any of the two series i F\ in (33) must reduce to a polynomial. 


Thus, depending on the value of d, three different cases are identified: 

(i) For — | < 9 < | the energy levels of the complex oscillator turn out to be given by 


K(e) = [n + - 


je 


n = 0,1,2, ••• , 


and the corresponding eigenfunctions read 


<j> n ( x ) = C„H n (^wx)e 2- 


w — e 2 , 


(34) 


(35) 


where C n are normalization factors and H n (i/u;a;) are Hermite polynomials of the 
complex argument y/cox. It is worth to note that the Hermite polynomials of complex 
argument appear in several physically interesting situations, e.g., in the calculation of 
the photon distribution for squeezed states and their Wigner distribution function 
in the eigenfunctions analysis for an optical beam in the paraxial approximation 


52 


53 


among others 54 . Let us note also that the Hermite polynomials of complex argument, 


in general, are not longer orthogonal 53 


(ii) For | < 9 < it turns out that 

E n {6)= (n + 0 e i ( 0_7r ), 
while the corresponding eigenfunctions become 


n = 0,1,2, ■ ■ - , 


(X 


= D n H„ 


ux)e 


\ux 2 


-u) = e 


(36) 


(37) 


(Hi) For 9 — ±| there are not square-integrable eigenfunctions of H since now the potential 


corresponds to the repulsive oscillator 54 


Thus, the eigenvalues of the complex oscillator potential lie always either in the first or in 
the fourth quadrant on the complex energy plane. This result is consistent with the fact 
that our potential is invariant under the symmetry transformation u j —y —u and, therefore, 
the eigenfunctions and eigenvalues for u ; and —u have to be the same, which is indeed the 
case. As a consequence, we can restrict ourselves to the domain — | < 6 < |. Moreover, 
since 

E n \-9) = [E n {e)}\ (38) 

it is enough to work with the domain 6 G [0, f), which we will do from now on. 

In analogy to the harmonic oscillator, it is possible to introduce now annihilation and 
creation operators a(j, which define a fixed direction u ; on the complex energy plane, in the 
following way: 

± 1 ( d 

a " := 7f( T d^ + ^ 


( 39 ) 











Since u> G C, then a+ and a u are not in general adjoint to each other. The algebra of these 


operators is analogous to the harmonic oscillator one, i.e., the Heisenberg-Weyl’s algebra 55 


[a“, a+] 

= u, 

(40a) 


= 2H 0 , 

(40b) 

[#o, a±] 

= ±wa^. 

(40c) 


Let us note that the eigenfunctions and eigenvalues for the complex oscillator can be found 
in an elegant alternative way by using the operators a±, as for the harmonic oscillator. The 
square-integrable eigenfunctions turn out to be 


</>n (%) 


-7= (a) 

Vn! 

C n H n (^x)e-> 2 , 


(41) 


where </>o oc e 


1 9 

~h ux 


is the ground state. Note that when e 


solution (33) can be reduced to (41), i.e., u (. x , E n ) = cf) n (x). 


E n ( 9 ) =: E n , the general 


A diagram of the complex energy plane, the eigenvalues E n of the complex oscillator for 
a fixed 6 as well as the action of the annihilation and creation operators afh can be seen in 
Fig. § 



Re ( E ) 


Figure 2: Diagram of the complex energy plane for a fixed 6 . Some eigenvalues are shown 
at their right positions. 

Let us apply now a non-singular fc-th order supersymmetric transformation which creates 
k new levels additional to E n for the new Hamiltonian HR = — + 14, i.e. 

Sp (H fc ) = {ej, E n ] j = 1, • • • , k) n = 0,1, 2, • • • } . 

In order to do that, we need to take either k solutions {/3i (x, £j) ; j = 1 ,k} of the 

initial Riccati equation for k different complex factorization energies Ej, or k solutions 
{Uj ; j = 1, • • • , k} of the corresponding Schrodinger equation. Since the oscillation the¬ 
orem is not longer valid for solutions of the stationary Schrodinger equation associated to 
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complex potentials, then the factorization energies {eq, • • • , £*,} can be chosen essentially at 
any position on the complex energy plane. This fact will be quite useful later on, when 
we will explore some particularly interesting (and simple) SUSY transformations. The new 


potential 14 is given by Eq. (16) with V 0 (x) = \lo 2 x 


The algebra for systems described by hh can be built departing from the complex oscil¬ 
lator one. Indeed, the natural ladder operators associated to hh are given by 


^2k+\ — > 


(42) 


where the fc-th order differential intertwining operators B * are given in Eq. (15). It turns 
out that the standard commutation relations of the PHA are satisfied: 


[Hfc,£ 2fc+1 ] — ±cu£ 2fc+1 , 


(43) 


which differ from Eq. (20a) in the factor u. In order to stick to the convention of Eq. (20a), 


let us redefine the Hamiltonian as 


H k 

TT _ K 

n-k — —, 


so that the relationship (43) is transformed into 


[H t , +1 ] = ±£± 


2fe+l' 


(44) 


(45) 


The Hamiltonians Hk associated with this algebra take the form 


H, = 


1 d 2 


+ —14 (x) = 


l d 2 


+ ^X 


1 d 2 


2 cj dx 2 l o 2u da: 2 2 cu dx 2 

which acquire the standard Schrodinger form through the substitution 

y = y/ujx. 


[In W (ui, • • • , Ukj ], (46) 


(47) 


The (2k + l)-th order ladder operators £^, +] and the Hamiltonian Hk (in fact a family) 
generate a 2/c-th order PHA. As we saw in section [3j systems connected with the PIV equa¬ 
tion must have third-order ladder operators and satisfy a second-order PHA. Thus, we need 
to identify the subfamily of Hk (if any) suitable for the purposes of this paper, i.e., having 
as well third-order differential ladder operators. In fact, there is a reduction theorem that 
will be useful later on, enabling us to reduce to second- the initial order (2k) of the PHA 


and thus leading to solutions of the PIV equation 37 . 


Let us analyze next some novelties arising in the complex case as compared with the real 
one. We will start by using the simplest situation, with k — 1, for illustrating them. Let 
us note that the transformation function associated to the factorization energy eq will be 
denoted as u\(x) throughout this Section, It will be taken either as one of the bound state 


eigenfunctions of Eqs. (35), (37) or as the general solution of Eq. (33). 
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4.1 First-order transformation 


As a first atypical fact, let us use as transformation function a bound state eigenfunction 
fa (x) of the complex oscillator. Note that the parity of fa is (—1)*, and the even eigenfunc¬ 
tions (with i even) do not have zeros on the real axis so they induce non-singular first-order 
SUSY transformations while the odd ones (with i odd) have one node at x — 0, inducing 
transformations with a singularity at this point. Let us study next each of these two cases 
separately. 


4.1.1 Bound state even seed solution 

Let us take as transformation function an even eigenfunction of H 0 : 

U\ (x) = H 2 j(y/ux)e~^ ux2 , j = 0, 1 , 2, • • • . 

The first-order SUSY partner potential of the complex oscillator turns out to be 

T r / \ 1 2 2, O' in- -n H 2 j —2 (\/Ukc) ' 2 

Vi (x) = -u x +u}~ 8 juj ( 2 j - 1 ) J - 2j J 

2 Yi 2 j(y/ux) V H 2 j (y/ujx) 

In Fig. [3] it is shown the typical behavior of V\ for several values of j. 


(48a) 


(48b) 




(a) 


x 

(b) 


Figure 3: Real (a) and imaginary parts (b) for the first-order SUSY partner potential V\ 
generated by an even eigenfunction of H 0 . Several non-trivial cases are shown for j = 1 
(blue), j = 2 (red) and j = 3 (dashed) with 9 — |. 


The explicit expressions for the eigenfunctions of Hi are given by 
,H 2 j_i (y/ux) 


{X) OC 


4.7- 


H 2 j(y/ux) 
so that the spectrum of H i becomes: 


H n (\4Jx) — 2nH n _i(\/a;x) 


e“5 wx , n = 0,1,2, ■■ 


Sp(Ri) = <J ( n + - ) u] n = 0,1, 2, • • • ; n ± 2 j J> , 


; n ± 2 . 7 , 

(49) 

(50) 
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i.e., the eigenvalue E 2 j of H 0 has been deleted in order to generate Hi, since the wavefunction 
annihilated by A[ in Eq. (10) is not square-integrable. 


As can be seen, for j > 0 the spectrum of Hi is composed of two equally spaced energy 
ladders: a finite one of length 2 j, which starts from Eq = ^ and ends at U 2 j_i = (2 j — \)“i 
an endless ladder departing from E 2 j+i = (2 j + |) u. Let us note that the lower ends of 
both ladders supply us with two extremal energies as well as the associated extremal states 
for our system. 

In Fig. [4] it is compared the normalized probability densities for the eigenfunctions of H 0 
with the corresponding transformed ones for H\. The graphs correspond to the case when 
the transformation function is the first even excited state (j) 2 ( x ) with 9 — |. 



(a) 

Figure 4: Comparison between the normalized 
of H 0 with n — 0 (blue), n = 1 (red) and n 
n — 0 (blue), n — 1 (red) and n = 3 (dashed), 
excited state (ft 2 ( x ) with 6 = |. 



x 

(b) 


probability densities for: (a) the eigenstates 
= 2 (dashed); (b) the transformed ones for 
The transformation function is the first even 


4.1.2 Bound state odd seed solution 

Suppose now that the transformation function is an odd eigenfunction of Hq, 

U\ ( x) = H 2 j+i( \/ujx ) e~^ wx2 , j = 0,1, 2, • • • , (51a) 

which vanishes at x = 0 and thus induces a singularity at the same point. Despite this fact, 
let us use the standard formulas of the first-order SUSY QM and analyze then the obtained 
results. 


The first-order SUSY partner potential reads now 


Vi (x) = 2 Uj2x ' 2 + w — 4 (2j + 1) oj 


2 j 


.H 2j -i(v / wic) 


H 


2j+l 


UJX ) 


(2j + 1) 


R 2j (y/ux) 

H 2 j+l(Vvx) 


(51b) 
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while the possible eigenfunctions of H\ are given by 




(i) 


X) oc 


(4 j + 2 ) H 2 j ~R n (y/ujx) - 2 nH n _i( v / wx) 

H 2j+1 (vW) 


e -> 2 


n = 0,1,2, ••• ; n ^ 2j+l. 

(51c) 


The potential V\ ( x ) has a singularity at x = 0, which can be interpreted as coining from 
a complex oscillator potential with an infinite barrier at the origin. Hence, the coordinate 
domain has to be redefined between two equivalent options: either (—oo,0) or (0, oo). For 
convenience, from now on we will take the domain x G (0, oo). 

The complex oscillator potential with an infinite barrier at the origin has odd parity eigen¬ 
functions, which vanish at x = 0 and x = oo. On the other hand, the transformed eigen¬ 
functions (51c) do not always satisfy the same conditions [56], thus we need to select the 


appropriate subset. In fact, it turns out that the yy ^ with n odd vanishes at x = 0 and 
x = oo while for n even it diverges at x — 0. Thus, the eigenfunctions of H\ in the restricted 
domain ( 0 , oo) are indeed: 


] (z) OC 


(4 j + 2 ) - (4 n + 2 ) H 2 n (\/wi) 


H 2i+ i(vW) 

The spectrum of Hi becomes now: 


e s"* , n = 0 , 1 , 2 , ••• ; n ^ j. 

(52a) 


Sp(tfi) = \E n = (2 n+-)u; n = 0,1, 2, • • • ; n ^ j } , 


(52b) 


i.e., the eigenvalue ( 2 j + oo of the complex oscillator with an infinite barrier at the origin 
has been deleted for generating the potential Vi (x), which domain is now x G ( 0 , oo). 


Similarly as for the non-singular case of section |4.1.1[ the spectrum of Hi can be seen 
as composed of two equally spaced energy ladders, but the spacing is now twice the original 
one for the non-singular case. In addition, the finite ladder has just j steps, starting from 
Eq = | oj and ending at Ej_i = ( 2 j — 4 ) oj. On the other hand, the infinite ladder starts 
from E j+ i = (2 j + |) cu. 

In Fig. [5] the typical singular behavior of Vi for several values of j is shown while in Fig. 
[ 6 ] the normalized probability densities for the eigenfunctions of H 0 and H\ in the restricted 
domain x G (0, oo) are compared. 


4.1.3 General seed solution 


For a transformation function U\ (x) = u(x,£i) of the form given in Eq. (33) with £i G C, 


the first-order SUSY partner potential of the complex oscillator turns out to be: 

2 


Vi (x) = ^ 2 x 2 


u'i (x) 


+ 


Ui (x) [ Ui (x) 


U 1 \?) 


(53) 


For some particular values of £i there is a further simplification since the confluent hyper¬ 
geometric function reduces to a polynomial. The simplest example appears for sq = — j for 
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(a) * (b) ' 

Figure 5: Real (a) and imaginary parts (b) for the singular first-order SUSY partner potential 
V\ generated by an odd eigenfunction of H 0 . Several cases are shown for j — 1 (blue), j — 2 
(red) and j = 3 (dashed) with 9 — 




(a) (b) 

Figure 6: Comparison between the normalized probability densities for: (a) the eigenstates 
of H 0 with n = 0 (blue), n — 1 (red) and n = 2 (dashed); (b) the transformed ones for n = 0 
(blue), n = 2 (red) and n — 3 (dashed). The transformation function is the odd excited 
state 03 (x) with 9 = 


which 


m (x) = e ^ x2 


In this case the new potentials become: 

1 


Ui (x) = ^ uj 2 x 2 


UJ 


1 H—^erf (y/cux) 

CO 


dx 


2 m 


(54) 


(55) 


ysa + ^erf(v^)i)J ’ 

which are complex extensions of the Abraham-Moses-Mielnik potentials derived several years 
The corresponding spectrum is composed of the levels E n = (n + |) u>, but now 


ago 

41 

57 

n = 


,0, 


n — — 1, 0,1, • • •. A plot of the real and imaginary parts of V\ for several values of are 
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shown in Fig. [7J The associated normalized probability densities for some eigenstates with 
a fixed value of zq are plotted in Fig. [8j 




(a) (b) 

Figure 7: Real (a) and imaginary parts (b) for the potential V 1 of Eq. (J55|. Several cases 
are shown for zq = — 0.6 + 0. 3i (blue), zq = 0.3i (red) and zq = 0.6 + 0.3z (dashed). We have 
taken £\ = —| and 6 — |. 




(a) * (b) ^ 

Figure 8: Comparison between the normalized probability densities for: (a) the eigenstates 
of H 0 with n — 0 (blue), n — 1 (red) and n = 2 (dashed) and (b) their corresponding 
quantities for H\. We have taken £i = —j, 6 = | and zq = 0.6 + 0.3i 


4.2 Higher-order SUSY partners 

A fc-th order SUSY transformation requires to choose in general, k independent seed solutions 
of the Schrodinger equation. In order to induce the reduction theorem for linking with the 
PIV equation, however, we will consider just connected seed solutions satisfying: 

uj (. x ) = (aj) J 1 u 1 (x ), 

U = - U ~ 3 = h"-,k, 
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(56a) 

(56b) 










with u\{x) being a complex seed solution associated to £\. As in the hrst order case, u\(x) 


can be taken either as one of the eigenfunctions of Eqs. (35), (37), as the general Schrodinger 


solution of Eq. (33), or as the solution given in Eq. (54). In particular, for k = 2 the second- 


order SUSY partner potential of the complex oscillator, for a fixed factorization energy £ 1; 
acquires the form: 


-r r / \ 1 2 2 W" (W' 

~W + \w 


(57) 


where 


W W(iii, a u Ui) oc ( oj 2 x 2 + u — 2 e ± ) u\ — ( u 'J 2 . (58) 

In Fig. [9] we have plotted the real and imaginary parts of V 2 , for the seed solution of Eq. ( [54] ) 
with £1 = — | and several values of zq. Some eigenfunctions of H 2 , {^n \ n = —2, —1, 0, • • • }, 


for £1 = — and zq = 0.9 + 0.4i, are shown in Fig. 10 




(a) 


(bf 


Figure 9: SUSY partners V 2 of Eq. (57). Real (a) and imaginary parts (b) for several values 
of the parameter zq: zq = 0.1 + 0.4i (blue), = 0.5 + 0.4f (red) and zq = 0.9 + 0.4z (dashed). 
We have taken £1 = — | and 6 — |. 


5 Connection with the Painleve IV equation 


We have seen previously that a second-order PHA with third-order ladder operators describes 
systems which are linked with the PIV equation (see section [3]). Thus, a first-order SUSY 
transformation applied to the complex oscillator provides solutions to the PIV equation. In 
fact, from Eq. (30c) we can obtain solutions to such equation in terms of the extremal states 
of the system by making carefully the changes from the harmonic to the complex oscillator: 





[In tp £j (a;)]', 


J = 1,2,3, 


(59) 


where the subscript is used to indicate the order of the SUSY transformation. The labeling 
of the wavefunctions in Eqs. (30) is not essential, so that by making cyclic permutations 
of the indices we can generate three different solutions to the PIV equation. Do not forget 
that U\ is a solution (not necessarily physical) of the Schrodinger equation (j7|) associated to 
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(a) 


(bf 


Figure 10: (a) Real part of V 2 (dashed) together with the real part of ifin'* for n = 0 (blue), 
n = 1 (red) and n — 2 (green). The magenta and purple lines represent the real parts of 
the new eigenstates V’-i (= ipe^) and (= V^) respectively, (b) The same as in (a) but 
for the imaginary part of V 2 and the corresponding eigenfunctions. We have taken E\ = — 

6 = | and ui = 0.9 + 0.4i. 


the factorization energy £\. Thus, for hxed value of £\ the solutions to the PIV equation 
constitute in fact a 3-parametric family. The three extremal states (up to constant factors) 
and their corresponding energies are given by 


Vfci oc Af^ 0) oc Afe ^ x2 , £i = 

Tp£ 2 OC —, S 2 = £l, 

Ml 

V ’£3 °C A l a u u 1) £3 = £1 + W. 


(60a) 

(60b) 

(60c) 


In Fig. 11 we show a set of complex solutions to the PIV equation for several values of 
the factorization energy and a hxed v. 




(a) 


(b) 


(2) (o') 

Figure 11: Solution g\ ' to the PIV equation. Real (a) and imaginary parts (b) of g{ ; for 
£1 = 0.01 + i (blue), £i = l + i (red) and £1 =2 + i (dashed). We have taken 6 — | and 
v = 0.8 + 0.5i 
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If "0£v ^ 0 V x € R, the functions g[ l] are non-singular. Furthermore, they have a null 
asymptotic behavior, g[ J> —> 0 when |x| —>■ oo. This can be seen from a different viewpoint in 
Fig. 12 for large values of |x|, the curve tends slowly to close around Re(gj^) = Im(g[ 2 ' ) ) = 0. 



Re(Si) 


(o') (o') (‘2') 

Figure 12: Parametric plot of Re(gq ) versus Im(c/) '). The function g\ ; tends to (0,0) in 
the limits |x| —> oo. We have taken e = 2 + i, 6 * = |, u = 0.8 + 0.5z and |x| < 10. 


For higher orders of SUSY (A: > 1), the solution to the PIV equation appears after the 
reduction theorem is produced. In fact, one of these solutions is given by 


9 k {x) 



1 W («!,■■■ ,ttfc-i) 

W (u 1 , • • • ,u k ) 


k> 2, 


(61) 


where the Uj, j = 1, • • • , k satisfy Eqs. (56). In Fig. 13 we illustrate the typical behavior of 
<72 for several values of eq and V\ = 0.8 + 0.5i 


It is worth to note that in the harmonic oscillator limit, Im(ei) = Im(z/i) = 9 = 0, it is 
possible to recover the results reported in |37|, where the number of nodes of the non-singular 
g functions for £i < \ are fixed by the order of the SUSY transformation, and they tend to 
zero when lx I —y oo. 


6 Conclusions 

In this paper we have introduced a complex generalization of the standard factorization 
method to analyze the SUSY partners of the non-Hermitian Hamiltonian associated to the 
complex oscillator. The corresponding potentials were generated through SUSY transfor¬ 
mations involving complex factorization energies. Some interesting transformations, which 
cannot be implemented in the real case without introducing undesired singularities in the new 
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(a) 


(b) 


Figure 13: Solution g 2 to the PIV equation. Real (a) and imaginary (b) parts of g 2 for 
£\ = 0.01 + i (blue), e 1 — 1 + i (red) and eq =2 + i (dashed). We have taken 6 — | and 
u\ = 0.8 + 0.5i. 


potentials, were found. In particular, the first-order SUSY transformation that uses bound 
state eigenfunctions can be implemented to obtain either non-singular potentials whose en¬ 
ergy levels become the initial ones (but for the employed bound state energy) or potentials 
with a singularity at x = 0, which are almost isospectral to the complex oscillator with 
an infinite barrier at the origin. In addition, it was identified a complex extension for the 
Abraham-Moses-Mielnik potentials, which turn out to be isospectral to the initial complex 
oscillator (up to an energy displacement). 

Let us note that the physical interpretation of complex energies remains open, although 
possible applications related to absorbent (dissipative) systems have been noticed recently 
0 (see also [53]). 

In order to finish, several additional open problems can be pointed out. In the first place, 
it is not clear the way we have to move on the complex energy plane in order to obtain a 
natural classification of the complex SUSY transformations, the associated potentials and 
the corresponding solutions to the PIV equation. In addition, this time we did not obtain 
any heuristic criterion to link the number of zeros of the PIV solution with the order of 
the SUSY transformation. We are convinced that these two subjects, although difficult, are 
worth to be explored in the future for clarifying better the peculiarities of the complex SUSY 
transformations discussed in this paper. 
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